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ABSTRACT 

We present new high-angular resolution observations toward the driving source of 
the HH 80-81 jet (IRAS 18162-2048). Continuum emission was observed with the Very 
Large Array at 7 mm and 1.3 cm, and with the Submillimeter Array at 860 fj,m, with 
angular resolutions of ~ O'.'l and ~ 0/! 8 respectively. Submillimeter observations of the 
sulfur oxide (SO) molecule are reported as well. At 1.3 cm the emission traces the well- 
known radio jet, while at 7 mm the continuum morphology is quadrupolar and seems 
to be produced by a combination of free-free and dust emission. An elongated structure 
perpendicular to the jet remains in the 7 mm image after subtraction of the free-free 
contribution. This structure is interpreted as a compact accretion disk of ~ 200 AU 
radius. Our interpretation is favored by the presence of rotation in our SO observations 
observed at larger scales. The observations presented here add to the small list of 
cases where the hundred- AU scale emission from a circumstellar disk around a massive 
protostar has been resolved. 
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1. Introduction 

It is well known that the formation of solar-type stars takes place with the assistance of an 
accretion disk that transports gas and dust from the envelope of the system to the protostar, and a 
jet that removes angular momentum from the system, allowing accretion to proceed (e.g., McKee 
& Ostriker 2007). While it is tempting to think that these mechanisms work all the way up in the 
stellar-mass range, it is not clear to what extent this assertion may be correct. Feedback from a 
growing protostar increases rapidly with mass, and early calculations suggested that it could be 
fatally disruptive for stellar masses M± > 8 M© (Wolfire &; Cassinelli 1987). It has been proposed 
that stars of higher mass are able to form via accretion when revised dust opacities and high mass- 
accretion rates are considered (Osorio et al. 1999). However, it is still unclear if processes such as 
radiation pressure, ionizing radiation, and jet/outflow feedback can terminate accretion onto the 
most massive young stellar objects, and if they do, how and when does this happen (for a review 
see, e.g., Zinnecker &; Yorke 2007). 

Theory suggests that several mechanisms can be at work to aid accretion onto massive young 
stellar objects (MYSOs). Flattened accretion flows focus the disruptive effects of radiation pressure 
(e.g., Yorke &: Sonnhalter 2002, Kuiper et al. 2011), photoionization (e.g., Peters et al. 2010), and 
jets and outflows (e.g., Wang et al. 2010, Cunningham et al. 2011) to some preferential angles while 
permitting accretion from other directions. Observations support this view; there are several known 
cases of flattened, dense gas structures that appear to be rotating and infalling around MYSOs 
(e.g., Cesaroni et al. 2007). However, these structures usually do not look like the Keplerian, stable 
disks with sizes of <100 AU seen around solar-type YSOs. Instead, they often appear very large, 
up to ~ 10 4 AU in diameter, and very massive compared to their central star(s), which therefore 
renders them unstable to fragmentation. In addition, they show apparent infall motions that are 
comparable in magnitude to their rotation. Indeed, model fitting of the SEDs of a sample of massive 
protostar candidates (De Buizer et al. 2005) suggests that the observed large-scale structures with 
sizes of thousands of AU could be naturally explained as infalling flattened envelopes, while the 
formation of the " true" accretion disks is expected to occur at scales of the order of the centrifugal 
radius (a few hundred of AU for these MYSOs). 

A more direct approach to test the presence of circumstellar disks in MYSOs is to look for 
"clean" examples of systems composed of a compact disk and a jet. These cases appear to be 
quite rare. Examples are G192. 16-3.82 (Shepherd et al. 2001), AFGL 2591 (Trinidad et al. 2003), 
IRAS 18162-2048 (Gomez et al. 2003), Cepheus A HW2 (Patel et al. 2005), IRAS 20126+4104 
(Hofner et al. 2007), IRAS 16547-4247 (Franco-Hernandez et al. 2009), and IRAS 13481-6124 
(Kraus et al. 2010). However, up to now, only Cepheus A HW2 has been observed with enough 
angular resolution and sensitivity to angularly resolve the emission of the disk; in the other cases, 
the disk emission has not been well resolved, or evidence for the accompanying jet is weak. 

One of the best candidates to look for a circumstellar disk is the massive protostar IRAS 
18162-2048. This protostar has a bolometric luminosity of L ~2xl0 4 L (Aspin & Geballe 1992), 
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equivalent to that of a BO zero-age main sequence star (M* > 10 M ). It powers a highly-collimated 
radio jet that extends 5.3 pc (at an adopted distance of 1.7 kpc; Rodriguez et al. 1980) toward 
the Herbig-Haro objects HH 80-81-81N (Marti et al. 1993, 1998). The jet is surrounded by a 
bipolar cavity seen at 8 /im (Qiu et al. 2008). Furthermore, it has been found that this jet is being 
collimated by a large-scale helical magnetic field, most probably originated in a rotating accretion 
disk (Carrasco-Gonzalez et al. 2010). Gomez et al. (2003) reported unresolved observations of 
the mm thermal dust emission from the exciting source of the jet. Recently, Fernandez-Lopez et 
al. (2011a,b) presented (sub)mm observations down to an angular resolution ~ O'.'S (~850 AU), 
and interpreted the emission as arising from a compact (size < 600 AU) accretion disk orbiting a 
~ 15 Mq central source. 

In this Letter, we present new sensitive observations performed with the Submillimeter Array 
(SMA) and the Very Large Array (VLA) towards the IRAS 18162-2048 MYSO. These observations 
resolve, for the first time, the dust emission of this utmost important source at angular resolutions 
down to ~ OVl, equivalent to ~ 170 AU. 



Observations in the 0.8-mm band were performed with the SMA[j (Ho et al. 2004) during two 
runs (2006 June 13 and 22). In the first run, the array was in its extended configuration, while in 
the second run, the array was in the compact configuration. Two sidebands spanning the frequency 
ranges 342.6-344.6 GHz and 352.6-354.6 GHz were covered. Calibration was performed using the 
MIR data calibration program. Quasars 3C454.3 and J1924-292 served as bandpass and phase 
calibrators, respectively. The absolute flux scale was derived from observations of Callisto, and is 
accurate to better than ~ 15%. Further processing and imaging was done in MIRIAD, AIPS, and 
IDL. In addition to the continuum emission, we report on the detection of compact emission of 
SO 8(8)-7(7) (^0=344.31061 GHz) around the core of the radio jet. HCN (4-3) and CS (7-6) were 
also detected but their emission is not confined to the immediate surrounding of the exciting source 
and their analysis is out of the scope of this Letter. 

The final continuum map was done from the extended configuration data with a uniform 
weighting to maximize the angular resolution (synthesized beam = O / .'97x0 / . / 70, P.A.= —26°). 
The SO map was made from the concatenated compact+extended data with an intermediate 
weighting (robust=0) to have the best compromise between resolution and sensitivity (synthesized 



2. 



Observations 



2.1. 



SMA Observations 



beam = l / /20xO / . / 99, P. A. 



34°). 



1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. 
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2.2. VLA Observations 

Observations at 1.3 cm and 7 mm continuum were carried out using the VLA of the National 
Radio Astronomy Observatory (NRAO^Jin its A configuration on 2004 November 7 (1.3 cm and 
7 mm) and 20 (1.3 cm). Phase and flux calibrators were J1820— 254 and 3C286, respectively. 
Data editing and calibration were performed using the AIPS package, following the standard high- 
frequency VLA procedures. 

Maps at 1.3 cm and 7 mm were made applying a tapering of 1750 kA and 2100 kA, respectively, 
in order to emphasize extended emission. Synthesized beams are / .'19x0 / . / 13 with a position angle 
(P.A.) of 20° (1.3 cm) and / .'12x0'. / 09 with a P.A. of 26° (7 mm). 

In order to compare the SMA observations with the emission of the HH 80-81 radio jet at 
similar scales, we calibrated VLA A configuration archive data at 3.6 cm continuum obtained in 
5 epochs (1990.2, 1994.3, 1995.5, 1997.1 and 2006.4). The map shown in this paper was made by 
concatenating data from all the epochs (synthesized beam=0 / .'5xO'. / 3; P.A.=0°). 



3. Results and Discussion 

3.1. An infalling rotating molecular envelope 

In Figure [T£i we show a superposition of the VLA map at 3.6 cm (contours) over the first 
moment of the SO molecule emission (colors) obtained with the SMA. At 3.6 cm, the radio jet 
appears with an elongated morphology along a P.A. of 20°. The SO molecule emission shows 
an extended envelope (size~3000 AU) around the driving source of the radio jet, with a velocity 
gradient roughly perpendicular to it (see Fig. [iji), that we interpret as rotational motions. In 
Figure [TJa, we show a position-velocity diagram along a direction perpendicular to the jet. From 
this diagram we measure a velocity gradient of ~ 2.5 km s _1 arcsec -1 , from which we infer a 
rotation velocity ~2 km s _1 at a radius ~ 1500 AU, assuming an inclination angle of 90° (i.e., the 
HH 80-81 jet is almost in the plane of the sky). The centrifugal radius (i.e., the largest radius on the 
equatorial plane that receives the infalling material) is given by R c = TqVq/ \GM), where vq is the 
rotation velocity at a distant reference radius ro, and M is the central mas^J Adopting the values 
of ro and i>o derived from our SO observations, and assuming M ~ 15 Mq (Fernandez-Lopez et al. 
2011b), we obtain R c ~ 650 AU. This result strongly suggests that the SO emission, which arises 
from radii larger than R c , is tracing an infalling and rotating envelope, while the "true" accretion 
disk should be formed at smaller scales, within the centrifugal radius. 



2 The NRAO is a facility of the National Science Foundation operated under cooperative agreement by Associated 
Universities, Inc. 

3 This expression for R c is derived assuming conservation of the specific angular momentum. Angular momentum 
losses during the infall process would decrease the actual value of the centrifugal radius. 
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3.2. A compact dusty disk 

In Figure^ we show a superposition of the 860 fim continuum map (contours) over the 3.6 cm 
map of the radio jet (colors). The compact 860 /im continuum emission is observed towards the 
core of the radio jet and has a flux density of 580±10 mJy. A Gaussian fit to the 860 /iin source 
gives a deconvolved FWHM <0!!7 (<1200 AU). Given that the spectral index a (where S u oc v a 
) of the free-free jet is ~ 0.2 at cm wavelengths (Marti et al. 1993), the free-free contribution at 
860 fim should be < 3 mJy (0.5 %). Therefore, the submm emission is dominated by dust, likely 
from an accretion disk (see below), but it remains unresolved at the SMA angular resolution of 
~0 / . / 8 (~1400 AU). 

In Figure]^ we show the superposition of our VLA 7 mm (contours) and 1.3 cm (color scale) 
maps covering the central region of the radio jet. The very high angular resolution of these maps 
(~0'. / 15, equivalent to 250 AU at 1.3 cm; ~0'. / 10, equivalent to 170 AU at 7 mm) allows us to resolve 
the structure of the core of the source. The 1.3 cm emission shows a jet-like morphology, consisting 
of a bright central source and two weaker sources to the NE and SW. The global orientation of the 
1.3 cm emission is similar to that of the larger scale radio jet detected with a higher signal-to- noise 
ratio at 3.6 cm. However, the central 1.3 cm source is oriented at PA.~10° at scales of ~0'. / l-0 / .'2, 
while at larger scales (~0 / . / 5) the jet is oriented at PA.~20°. This suggests a precession of the jet 
axis, as previously proposed by Marti et al. (1993) for the large-scale (10') jet. 

Emission at 7 mm is detected only towards the center of the radio jet and shows a quadrupolar 
morphology which can be described as the superposition of two overlapping, elongated sources. 
We have fitted the 7 mm source with two Gaussian ellipsoids. This fit is shown in Figure [3] and 
the obtained parameters are listed in Table [1} Both components are extended with perpendicular 
orientations (roughly N-S and E-W) and similar flux densities (~3 mJy; see Table [l]). 

Similar quadrupolar morphologies have been observed in other radio jets of low-mass stars at 
wavelengths where comparable contributions of free-free emission from the radio jet and thermal 
dust emission from a perpendicular disk are present (e.g., HH 111 at 7 mm: Rodriguez et al. 2008; 
HL Tau at 1.3 cm: Carrasco- Gonzalez et al. 2009). The total flux density at 7 mm (~6.3 mJy) and 
the flux density of the central 1.3 cm component (see Table [T]) imply a spectral index a ~2, much 
higher than the spectral index of the ionized jet, strongly suggesting the presence of thermal dust 
emission at 7 mm. Therefore, we interpret the morphology of the 7 mm emission as the result of a 
combination of free-free emission from an ionized jet (N-S component) and thermal dust emission 
from a perpendicular structure (E-W component). In Figure [2j: we show a superposition of the 
E-W component (after substraction of the free- free contribution) over the 1.3 cm emission of the 
radio jet. From the deconvolved major axis of the E-W component (see Table[T]), we derive a radius 
of ~200 AU, consistent with the radius (;$300 AU) of the unresolved 1.4 mm continuum source 
detected by Fernandez-Lopez et al. (2011a) with the SMA. The radius obtained from our 7 mm 
observations is smaller than the centrifugal radius inferred from our SO observations {R c ~ 650 
AU). Thus, although our observations lack direct kinematical information, we interpret the 7 mm 
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emission elongated in the E-W direction as tracing a dusty accretion disk. 

The inferred radius of the disk (~200 AU) is similar to that of the rotating disk around 
the massive protostar Cep A HW 2 (~300 AU; Patel et al. 2005), and fully consistent with the 
theoretical estimates for MYSOs (De Buizer et al. 2005). It is also of the same order (although 
somewhat larger) than typical values of disks around low-mass protostars (<100 AU, e.g., Andrews 
et al. 2009). 

A rough estimate of the disk mass can be obtained from the dust emission at 7 mm using the 
equation 
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where n v is the dust opacity per gram of dust+gas at frequency v, S u is the flux density, T d is the 
dust temperature, and D is the distance to the region. The main sources of uncertainty in the mass 
determination are and k v (sensitive to the distribution of grain sizes and to the presence of ice 
mantles on the grains). Since the upper- level energy of the observed SO transition is 87.45 K, then 
the temperature at 1500 AU scales should be >90 K. Assuming that the temperature varies with 
the radius as T a R~ 5 , we obtain a lower limit of >250 K for the temperature at 200 AU (the 
radius of the disk). For the 7 mm opacity we adopt a value of 0.001 cm 2 g _1 (assumes a gas-to-dust 
ratio of 100), appropriate for warm (T ~ 300 K) gas with a grain maximum size of 1 mm (D'Alessio 
et al. 2001). With these assumptions and using Eq. 1, we estimate Md; s k <4 M . This would 
imply a ratio of the disk to stellar mass M^^/M^ <0.3 which is towards the upper end of typical 
values found for low-mass protostars (M^k/M* ~0.001-0.1; e.g., Andrews et al. 2009, Ricci et al. 
2010). Our estimate for the disk mass is consistent with the value given by Fernandez-Lopez et al. 
(2011a). However, it should be noted that both estimates are very uncertain since they have been 
obtained from very simple assumptions. In order to obtain an accurate estimate of the mass of the 
disk it would be necessary to model the SED over a wide range of frequencies using an accretion 
disk model that takes into account effects such as the geometry, physical properties, and heating 
from the protostar or the dusty envelope in a physically self-consistent mode. 

The accretion rate of the disk onto the star cannot be easily calculated from our present data. 
However, a lower limit to the accretion rate can be estimated because this should be higher than the 
mass-loss rate in the jet. We assume a pure hydrogen jet with constant opening angle 6q, terminal 
velocity Vj ct , ionization fraction xo, and electron temperature T e = 10 4 K. We further assume that 
the jet axis is in the plane of the sky. Under these assumptions, and following Reynolds (1986), 
the mass loss rate in the jet is given by 
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The opening angle of the jet is estimated to be #o=2 arctan(# m i n /0 ma j), where 9 m \ n and # ma j 
are the deconvolved minor and major axes of the jet, respectively. Using the deconvolved sizes and 
flux of the jet component from Table 1, taking Uj ct =1000 km s _1 (Marti et al. 1995), and assuming 
xo=0.1 (e.g., Rodriguez et al. 1990, Shang et al. 2007), we estimate a mass-loss rate M out ~10~ 5 
M yr _1 . The accretion rate M acc of the disk onto the star is expected to be ~10 times larger 
than the mass-loss rate (e.g., Bontemps et al. 1996), or M acc — 10 -4 M & yr -1 . This value of 
the accretion rate is ~ 100 times larger than typical values for low-mass protostars (M acc < 10~ 6 
M Q yr _1 ; Evans et al. 2009). 



4. Conclusions 

We found that the IRAS 18162-2048 massive protostar, that was known to be associated with 
the highly collimated jet HH 80-81, is surrounded by a compact disk that is angularly resolved 
in our observations. The observed radius of the disk (~ 200 AU, similar to the radius of the 
Cep A HW2 disk) is in agreement with the values theoretically expected for massive protostars, 
and somewhat larger than typical values observed in low-mass protostars. The ratio between the 
mass of the disk and that of the central star (~0.3) is within the range (although towards the upper 
end) of typical values found for low-mass stars. The mass accretion rate of the disk onto the star 
(~10~ 4 M yr _1 ) seems to be much higher than in the low-mass case. These results reinforce the 
idea that the formation of high-mass stars is governed by physical mechanisms that are similar (but 
scaled-up) to those of low-mass star formation. 
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Table 1. Parameters of continuum sources 



Wavelength 


Position (J2000) a 


Flux Density 




Deconvolved 




RA DEC 


(mJy) 




Angular Size b 


1.3 cm 


18 19 12.094 -20 47 30.91 


1.7 ± 0.2 


0723 


x 0707 ; 10° ± 7° 


7 mm (N-S) c 


18 19 12.094 -20 47 30.92 


2.7 ± 0.3 


0712 


x 0702 ; 11° ± 7° 


7 mm (E-W) c 


18 19 12.096 -20 47 30.90 


3.6 ± 0.5 


0723 x 0713 ; 110° ± 10° 


860 /jcm 


18 19 12.084 -20 47 30.84 


580 ± 10 




<077 



a Units of right ascension are hours, minutes, and seconds and units of declination are 
degrees, arcminutes, and arcseconds. The absolute positional accuracy is estimated to be 
0705. 

b Major axis x minor axis; position angle of major axis. Uncertainty in major and minor 
axis is estimated to be 0702. 

c Emission at 7 mm shows a quadrupolar morphology (see text). Parameters are obtained 
by fitting the 7 mm emission with two Gaussian ellipsoids using the task JMFIT. As initial 
guesses for the fit we used for both components the same position (that of the peak of 
the total 7 mm emission) and the same peak intensity (half of the peak of the total 7 mm 
emission). For the initial P.A.s, we used 10° (as inferred from the 1.3 cm source) and 100°. 
Fits using different initial values converge to final values within the uncertainties. 
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Velocity (km/s) 
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Fig. 1. — (a) Superposition of the 3.6 cm VLA contour map over the SMA first-order moment (velocity) color 
map of the SO molecule emission. Contours are —3, 3, 6, 12, 25, 50, 100, and 200 times the rms of the VLA map, 
10 /xJy beam" 1 . Color scale ranges from 11.0 to 14.5 km s _1 . (b) Position- Velocity diagram at a P.A. of 110° 
(perpendicular to the radio jet) centered on the core of the radio jet. Contour levels are 0.4, 0.8, 1.0, 1.4, 1.8, 2.2 and 
2.6 Jy km s _1 beam -1 . Spatial and velocity resolutions are represented by the cross in the lower-left corner. 
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Fig. 2. — (a) Superposition of the 860 /im continuum SMA extended configuration map (contours) over the 3.6 
cm VLA map (colors). Contours are —3, 3, 4, 8, f6, and 32 times the rms of the SMA map, 10.7 mJy beam" 1 , (b) 
Superposition of the 7 mm VLA map over the 1.3 cm VLA map. Contour levels are 3, 4, 6, 8, 10, 14 and 18 times 
the rms of the map, 120 /xJy beam -1 , (c) Same as (b), after subtraction of the free-free contribution at 7 mm. 
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Fig. 3. — Decomposition of the quadrupolar 7 mm source in two perpendicular Gaussian sources. Panels (a) and 
(b) show the VLA 7 mm map and the fitted model, respectively. Panel (c) shows a comparison of the 7 mm data 
and the model. Panel (d) shows the decomposition of the 7 mm emission in two perpendicular components. The N-S 
(respectively, E-W) emission is obtained by substracting to the 7 mm data the E-W (respectively, N-S) model. Panel 
(e) shows the two Gaussian model sources whose parameters are given in Table [l] We note a small difference in the 
central positions of the components that, if real, could be due to density inhomogeneities in the jet and/or the disk. 
Panel (f) shows the residual after subtracting the model to the 7 mm data. In all panels, contours are —3, 3, 4, 5, 6, 
8, 10, 12, 16, and 20 times the rms noise of the 7 mm map, 120 /iJy beam -1 . 



